Neutral Solar Wind Generated by Lunar Exospheric Dust at the Terminator 
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Abstract. We calculate the flux of neutral solar wind observed on the lunar surface at the terminator due to 
solar wind protons penetrating exospheric dust grains with (1) radii greater than 0.1/im and (2) radii greater 



than 0.01/xm. For grains with radii larger than 0.1/im, the ratio of the neutral solar wind flux produced by 
exospheric dust to the incident ionized solar wind flux is estimated to be ^10~ 4 -10~ 3 for solar wind speeds in 
excess of 800 km/s, but much lower (<10 -5 ) at average to slow solar wind speeds. However, when the smaller 



00 

grain sizes are considered, this ratio is estimated to be >10 5 at all speeds, and at speeds in excess of 700 km/s 
reaches ~10 -3 . These neutral solar wind fluxes are easily measurable with current low energy neutral atom 

C/3 . 



instrumentation. Observations of neutral solar wind from the surface of the Moon would provide independent 
information on the distribution of very small dust grains in the lunar exosphere that would complement and 



constrain optical measurements at ultraviolet and visible wavelengths. 



1. Introduction 



of small dust grains [e.g. Wimmer-Schweingruber and Bochsler, 2003]. 



Although the vast bulk of solar wind hydrogen is ionized, a small fraction of the solar wind observed at 1 AU 

> 

^vj and beyond in the inner hcliosphere is neutral. The neutralization likely occurs either through charge exchange 
IT)'. . 

, with other neutral atoms [Collier et al., 2003]; for example, exospheric or interstellar neutrals, or by penetration 

(N 

o\ t 

Despite neutral solar wind having been discussed for many years [Fahr, 1968; Holzer, 1977] it has only 
recently been observed by the Low Energy Neutral Atom (LENA) imager at Earth [Collier et al., 2001] and by 
the Analyzer of Space Plasma and EneRgetic Atoms (ASPERA-3) at Mars [Brinkfeldt et al., 2006]. At Earth, 
neutral solar wind generated by the charge-exchange interaction of solar wind protons with the exosphere has 



provided remote observations of changes in the location of the magnetosphere and high latitude reconnection 
points [e.g. Taguchi et al., 2005, 2006; Collier et al., 2005]. The difficulty observing neutral solar wind arises from 
both technical challenges associated with detecting low energy neutral atoms and the potential high background 
levels due to stray light that results from looking in a direction close to the Sun. In LENA's case, these problems 
were mitigated by conversion surface technology a light-tight design with a hemispherical analyzer, and a time- 
of-flight unit to differentiate particles from photons. In ASPERA's case, the observations at Mars were taken 
during a period when Mars Express was entering eclipse after solar occultation but before Mars blocked the 
neutral solar wind signal. 

Future neutral solar wind observations may be made on ESA's Solar Orbiter mission which could fly the 
first dedicated neutral solar wind analyzer [J. Hsieh, private communication, 2007]. 
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Although the solar wind contains a neutral component arising from many different sources, we focus here on 
the locally-generated lunar neutral solar wind as observed by a hypothetical neutral atom instrument similar to 
those flown on previous missions and deployed on the lunar surface in much the same way as the Apollo Lunar 
Surface Experiments Package (ALSEP) and the proposed lunar sortie science missions. 

2. Lunar Dust 

Without exception, the astronauts who walked on the surface of the Moon confronted problems due to 
lunar dust: it adhered to clothing and equipment, it reduced visibility, and it caused difficulty breathing [e.g., 
see Stubbs ct al., 2007a and references therein]. In fact, "the invasive nature of lunar dust represents a more 
challenging engineering design issue, as well as a health issue for [lunar] settlers, than does radiation" [Schmitt, 
2006]. Furthermore, lunar dust may pose acute toxicity risks to astronauts [Liu et al., 2007; Park et al., 2006]. 

Evidence from Surveyor vidicon images of horizon glow [Rennilson and Criswell, 1974], as well as excess 
brightness in photographs of the solar corona taken by Apollo astronauts just inside the Moon's shadow [McCoy, 
1976] and Apollo astronaut observations of horizon glow and "streamers" immediately prior to orbital sunrise 
[McCoy and Criswell, 1974], all suggest that a substantial population of exospheric dust exists in the terminator 
region of the Moon. Rennilson and Criswell [1974] concluded that ~10 7 more particles per unit time must be 
ejected than can be achieved by meteoritic impacts in order to produce a dust cloud that could scatter enough 
sunlight to create the horizon glow observed by the Surveyor landers. 

Solar photoelectron emission on the dayside tends to drive the lunar surface positive, typically about +10 V, 
while solar wind electrons incident on the nightside drives the surface negative, typically about -100 V [Freeman 
and Ibrahim, 1975; Stubbs et al., 2007b, c], although occasionally regions of the lunar surface can charge up to a 
few kilovolts negative [Halckas et al., 2007]. It is believed that differential charging of the lunar surface in the 
terminator region can result in strong local electric fields that are able to eject charged dust into the exosphere 
[e.g., Criswell and De, 1977: De and Criswell, 1977: Borisov and Mall, 2006; Wang et al., 2007]. 

From the limited observations of exospheric dust phenomena, it appears that large grains (for our purpose) 
of ^5-6 /im radius are "levitated" to heights of ~ 10-30 cm [Rennilson and Criswell, 1974] while the smaller 
grains, ~0.1 /im, which are more relevant to neutral solar wind formation, are "lofted" to heights in excess of 
^100 km [McCoy, 1976]. It is estimated that the exospheric dust observed by the Apollo astronauts had a scale 
height of ^10-20 km [Zook and McCoy, 1991]. Ultraviolet horizon glow observations by Page and Carruthers 
[1978] on the lunar dayside could be interpreted as being due to dust with radii as small as 0.01 /um. Recently, 
the Apollo samples have revealed evidence for lunar dust grains as small as 0.01 /im in diameter [Taylor, 2007; 
Liu et al., 2007]. 

The electric field configuration near the terminator region appears to drive the electrostatic transport of 
charged dust [Farrell et al., 2007; Stubbs et al., 2006, 2007c]. The most compelling evidence for this came from the 
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Lunar Ejecta and Meteorites (LEAM) experiment deployed on the surface of the Moon by Apollo 17 astronauts 
[Berg et al., 1976]. LEAM was designed to detect very small fluxes of hypervelocity (>10 km/s) impacts from 
interplanetary and interstellar dust, but instead detected significantly larger fluxes of slower moving (~100 m/s) 
highly charged dust grains of lunar origin [Berg et al., 1976; Colwell et al., 2007]. The limited observations of 
the lunar "dust-plasma" environment suggest that it is most active at the terminator, therefore we focus our 
attention on dust-created neutral solar wind near this region. In addition, this is also the region where this 
environment will most likely pose the greatest hazard to human and robotic explorers [Stubbs et al., 2007a; 
Farrcll et al., 2008]. 

3. The Effective Dust Cross Section 

Although neutral solar wind also forms through interaction with neutral atoms, we will calculate here only 
the flux of neutral solar wind created from the penetration of exospheric dust. Further, we assume that lunar 
dust grains are spherical with radius a, although lunar dust can actually be quite irregular [Carrier et al., 1991; 
Liu et al., 2007]. However, there has been no direct characterization of exospheric dust. We note that in the 
following calculations the irregular shape of lunar dust grains will tend to increase the flux of neutral solar wind 
because a larger fraction of the dust grains of a given mass will be penetrable. 

We consider solar wind protons, which range in speed from about 200 km/s (209 cV) to 1200 km/s (7.5 keV) 
with a nominal value of about 437 km/s (1 keV). If a solar wind proton impacting a dust grain has a range A 
(typically 0.01 (im to 0.1 ^m) in the dust grain, then provided A > 2a the proton can penetrate the dust grain 
regardless of where it hits (i.e. independent of its impact parameter). Thus, for A > 2a, the effective dust cross 
section is ira 2 , the projected area of the grain. 

However, if A < 2a, then the proton can only penetrate the fraction of the dust grain's projected area around 
the edge of the grain where its path length does not exceed its range. If the proton's path through the grain 
has a minimum distance to the grain's center (i.e. impact parameter) of r, then it must penetrate a length of 
2\/a 2 — r 2 to make it through the dust grain, such that the requirement for dust penetration is 



A > 2^/a 2 - r 2 



or 



(1) 




(2) 



This means that only the cross-sectional area between r = \J a 2 — (A/2) and r = a contributes to the 
effective cross-section such that the effective cross-section, <r e ff(a) becomes: 




(3) 



Thus 



Veff{a) 




for A < 2a 
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cr e ff (a) = tto? for A > 2a. (4) 

Interestingly, for larger dust grains (A < 2a), the cross section is independent of grain radius. This is because 
although larger grains have a greater total cross-section, a smaller fraction of that total cross-section can be 
penetrated such that the two effects cancel. 



4. The Lunar Dust Concentration at the Terminator 

Murphy and Vondrak [1993] used dust radii and altitudes (scale heights) from Rennilson and Criswell [1974], 
McCoy [1976], and Zook and McCoy [1991] together with the assumption that the scale height is a power law 
function of grain radius, to estimate that 

Z ° = ^873' ( 5 ) 

where zq is the scale height and a is the dust radius, both in meters, and C — 2xl0~ 15 m 11 / 3 . 

Following the work of Murphy and Vondrak [1993], we assume the dust concentration for a given grain 
radius above the lunar surface is a power law with an index of unity, and that exospheric dust has an expo- 
nential distribution with height above the surface (z) at a given dust radius. However, the estimates of lunar 
dust distributions from scattered light observations are still controversial, and should be considered "order of 
magnitude" as should the predictions presented here. However, note that all our assumptions in the paper tend 
to be conservative such that they are more likely to underestimate the actual flux of lunar dust-generated neutral 
solar wind. 

Thus, we take the dust concentration per unit dust radius at the terminator to be 

77-0 (X^^^Z 

p(a, z) = — cxp{ — } , (6) 

where altitude z is in meters and n is a scaling factor with dimension m~ 3 . 

Using expression (6), we can calculate the dust column concentration at the terminator T, for grain radii 
greater than some a m ; n as 

8/3, 



r°° r°° nQ a s/3 z 
/ / — exp{ — ) dzda 



3no£ 

8„ 8 / 3 ' 
Sa min 

Murphy and Vondrak [1993] determined T and hq using the sky brightness observed during lunar "twilight" 
by the astrophotometer aboard Lunokhod-II [Severny et al., 1974], together with the estimates of exospheric 
dust concentrations from McCoy's [1976] model "0". The best fit was found for T = 1.4xl0 9 m~ 2 for dust radii 
from 0.1 fim to 6 (im, as reported by Stubbs et al. [2007d]. 

Equating this value of T to expression (7), using 10~ 7 m (0.1 fim) as a m ; n , and noting that for a maximum 
dust radius of 6xl0~ 6 m (6 /xm) the integration to infinity is a good approximation, we get 

no = 4.0xl0 5 m" 3 . (8) 
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Thus, the expression we will use for the exospheric dust distribution is a modified version of the Murphy 
and Vondrak [1993] model, 

. , 4xl0 5 , a s / 3 z -, 

p[a , z) = _ r ^__Y (9) 

where p is in units of m~ 4 and a and z are in units of m. 

It is interesting to use expression (9) to calculate the concentration versus altitude profile n(z) for dust 
grains larger than some radius a min at the terminator: 



/•OO 

n(z) = / p(a, z) da 

J Qmin 

f°° 1 a 8 / 3 z 
= 4X1 ° 5 L ■ a eXP ^2^^ a ' 



where n is in units of m 3 . By using the change of variable 



2x10 



-15 ' 



this integral may be recast as 



3 f°° 1 

n(z) = - • 4xl0 5 J a s/3 z — exp{— x} dx 



8/3 



where 



nn\ 

n—l 



(10) 



(11) 



= 1.5xl0 5 E 1 ( mi " ), (12) 
H 2xl0- 15; ' 



Ex(x) =r(0,a:) 

= / cxp{-t}t~ 1 dt 

J X 

= -1'-i»M-£^, (13) 



where 7 = 0.577 is Euler's constant [Abramowitz and Stegun, 1972]. 
Thus, for a m ; n = 0.1xl0~ 6 m, we get 

n(z) = 1.5xl0 5 Ei(l.lxl0" 4 • z). (14) 

Therefore, as long as z<§;(l.lxl0~ 4 m -1 ) -1 = 9.3 km, we can use the approximation 

Ei (a) « - 7 -ln{x}. (15) 

Thus, 

n(z) = 1.5xl0 5 [- 7 -ln{l.lxl0" 4 z}] 

= 1.5xl0 5 [-ln{e 7 }-ln{l.lxlO" 4 z}] 

= 1.5xl0 5 ln{^}. (16) 
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Figure 1 shows the exospheric dust concentration as a function of altitude at the lunar terminator from 
Stubbs et al. [2007d] along with the approximation given by equation (16). The approximation appears to be 
reasonable within a factor of two or so, at least up to 1 km altitude. 

5. "Visible" Light-Scattering Dust Model 

In this model, we assume that the smallest lunar dust grains are those that are effective at scattering visible 
light (a min =0.1 fim). We will also assume the composition of lunar grains is silicon mainly because data on 
particle range through silicon are readily available. However, this seems a reasonable approximation as even 
mafic materials have a very high silica content. The composition of lunar dust does change with dust radius: 
crystalline silica and plagioclase/mafic mineral ratios are greater in the a < 5 /tm fraction than in the coarser 
fractions for all lunar soils. Comminution (breaking into fragments), agglutination (melting), vertical mixing, 
and local lateral transport control the composition versus radius relationship [Hcikcn ct al., 1991]. 

Figure 2, based on the data in figure 4a of Dcmond ct al. [1980], shows the range in silicon in ^m as a 
function of hydrogen speed in km/s. It is not anticipated that there is any substantial difference in the data 
between hydrogen and proton penetration. The figure suggests that a good approximation for the range of 
protons in (im through lunar dust is 

A = 8.9xl0~ 8 • v 2 
= 0.017 ^rn/keV, (17) 

where v is the proton speed in km/s. 

It should be noted that equation (17) is an average projected range and that the actual range for any given 
proton can vary considerably. We use this expression as a reasonable estimate with some protons penetrating 
farther and some not as far. Also, this will depend on the shape characteristics of exospheric dust which are not 
well-known. This means that, considering only "visible" light-scattering dust, all reasonable solar wind velocities 
fall in the regime for A < 2a (see equation 4) such that 

VeffW = —A 

= 6.2xlO~ 2 V, (18) 

where a e ff is in m 2 and v is in km/s. 

Based on the data from Figure 1 of Kallcnbach et al. [1993], at these energies, >80% of the protons exiting 
thin carbon foils are neutral. Although it is not clear whether these data are directly applicable to this problem, 
the charge state distribution of protons through lunar dust is likely similar [Wimmer-Schweingruber and Bochslcr, 
2003]. We will assume because such a large fraction of protons exiting carbon foils are neutral that all protons 
that penetrate lunar dust grains emerge neutral. 
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Figure 3 shows the viewing geometry, which is tangential to the lunar surface at the terminator. Here x is 
the distance along the line-of-sight viewing into the solar wind at the terminator, R is the lunar radius, 1738 km, 
and z is the height above the lunar surface of the point x. Thus, 



(z + R) 2 = x 2 + R 2 , 



(19) 



or 



z = R v / l + x 2 /R 2 -R 

„2 

1 



R 



2~R' 



x' 

2R 2 



R 



(20) 



where we have expanded to first order under the assumption that x <C R. This expansion is not formally correct 
for very small grains which have a large scale height. However, equation (20) is still sufficient for our purposes 
(sec Section 6). 

We can calculate the ratio of the neutral solar wind flux $NSWi to solar wind flux $ sw , by integrating over 
all dust radii along the line-of-sight [Roelof, 1997; Roelof and Skinner, 2000]: 

ff. roo r oo „8/3. 



J NSW 



\ I — exp{ — ) • (Jeff {a) dxda. 

Ja^in JO a L 



(21) 



In expression (21), we have used the dust distribution function given by equation (9) and taken the effective cross- 
section to depend only on the proton range because of the relatively large radius of the dust grains considered 
here: 



Thus, 



NSW 



Veff{a) = -A 2 



tt , 2 [ x n Q [°° a 8 / 3 x 2 



Noting that integrating the normalized Gaussian gives 



the x-integral above becomes 



Consequently, 



,8/3 



CX P{-— -2R^ dX 



$nsw 7i" 3/2 A 2 



n 



3^ 3 / 2 A 2 n fCR 
16 V 2 



27TCT, 



nCR 



-4/3 



a- J/3 da 



4/ 3 



(22) 



(23) 



(24) 



(25) 



(26) 
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Using the values A = 8.9xl0~ 14 v 2 m from equation (17), n = 4.0xl0 5 m -3 from equation (8), C = 
2.0xl0~ 15 m 11 / 3 , R — 1738xl0 3 m, and a m ; n = 0.1xl0 _6 m for the smallest dust grain radius considered here 
(which the solar wind is taken to not fully penetrate), we get 

^ = 3.0xlO- lf V. (27) 

<S>sw 

The solid line in Figure 4 shows this ratio plotted as a function of solar wind speed. At high solar wind 
speeds, even considering only "large" dust grains, the ratio of neutral solar wind to solar wind flux comes close 
to 10 -3 . Note that in this expression, because the solar wind does not completely penetrate the grains for any 
reasonable speed, we never reach a regime in which the neutral solar wind flux is directly proportional to solar 
wind flux. Interestingly, this non-linear behavior also exists for magnetosheath charge exchange at the Earth 
[Collier ct al., 2005] and may very well be a general property of neutral solar wind production. 

For reference, the dotted line in Figure 4 shows the neutral solar wind count rate that an instrument like 
IMAGE/LENA [Moore et al., 2000] would observe at the terminator assuming a hydrogen efficiency of 2xl0~ 4 , 
an aperture size of 1 cm 2 , and an average solar wind flux of 3xl0 8 /cm 2 /s. 



6. "Ultraviolet" Scattering Dust Model 

In this model, we assume that the smallest lunar dust grains are those that are more effective at scattering 
ultraviolet light (a min = 0.01 /zm). It is believed that grains of this radius are easily lost since they are prone to 
tribocharging and attach to everything. Hence, when the Apollo samples were transferred between containers, it is 
highly likely that the smallest grains were preferentially lost, which would result in their being under-represented 
in any subsequent size distribution analysis. However, UV horizon glow seen by Page and Carruthers [1978] in 
the 1050A to 1600A (105 to 160 nm) range could be interpreted as due to dust as small as 0.01 /im. Furthermore, 
grains this small are present in the Apollo samples [Park et al., 2006]. 

In order to fully penetrate the smallest grains in the population, A > 2a: 

2a = 0.02 (Um < 8.9xlO"V, (28) 

or 

v > 474 km/s = v crit , (29) 

so that below about average solar wind speeds, the solar wind does not fully penetrate any of the lunar dust 
grains, and we can use the formalism for the "visible" dust model, but of course with the neutral solar wind 
fluxes scaled upwards by a factor of IO 4 / 3 = 21.5 because the minimum grain radius in the "UV" model is a 
factor of 10 smaller than the "visible" model (0.01 /im versus 0.1 ^m) described by equation (26). Thus, 

= 6.4xl0~ 15 v\ v < 474 km/s (30) 
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For solar wind velocities in excess of v cr i t , the solar wind fully penetrates lunar dust grains with radii smaller 
than a tr , so-called because it marks the "transition" radius from partial to full dust penetration, where 



a tr =4.5xlO"V, 



(31) 



from equation (28) with v again in km/s and a tr in /an. For dust grain radii smaller than a tr , the solar wind 
fully penetrates the grains and the effective cross section 



a eff(a) = na 2 , a < a tr 



(32) 



while for grain radii in excess of a tr , the effective cross section is 



°eff(a) = —A 2 



• 7.9xlCT z 'v 



274 



a > a tr 



(33) 



with v in km/s. This is, of course, just Tra 2 r because the cross section is independent of grain radius. 
So, for the flux ratio at solar wind speeds in excess of v crit , we simply start with equation (21), 



NSW 



= n 



/ / — exp{ —)■ a e ff{a) dxda 

ia mln JO a *~ 

\ cr e ff(a) ■ a~ 7/3 da, 

Ja min 



ir£R 



(34) 



but note that the cross section, cr e //(a), now depends on grain radius with the dependence beginning at grain 
radii smaller than a tr : 



$nsw / 7r£i? 
= n o 



= n 7r 3 / 2 < 



ra tI poo 

/ na 2 ■ a- 7 ^ da + -\ 2 ■ a~ 7 / 3 da 

J a m i n J at r 

/ a-V3 da+^- a- 7 /3 da 

-J a. m in J atr 



= no7r 3/2j^.3r a 2/3|»trj 



7T 



3 / 2 A 2 n fCR 



2 2 
3tt 3 / 2 n /Zfl 
~2 



2/3 _ 2/3 
a tr a min "+" 



4 

A 2 1 



[a- 4 / 3 1 



8 at /3 



(35) 



This expression can be re- written as 

$nsw 3tt 3 / 2 n fCR A 2 1 



2 8 4/3 



A 2 



[a 



4/3 2/3 



I ^ &min ^ 



4/3' 



(36) 



where the expression outside of the large square brackets to the right of the equals sign is expression (26) and 
represents the flux ratio achieved if the protons do not fully penetrate the grains. The expression inside the large 
square brackets, then, may be viewed as a "correction term" which takes into account that the cross section does 
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not increase with increasing solar wind speed once the solar wind fully penetrates a particular grain size. Noting 
that A = 2a tr , the term in the square brackets C(a tr ), may be expressed as 

4/3 / N 2 

j _ 2 f^mi, 

a t i J V a tr 



C(a tr )= 3(^1 -2(^1 , (37) 



^NSW 


= 6.4xl0" 15 w 4 


$sw 




^NSW 


= 6.4xl0~ 15 v 4 







which, because at r > a m in, C(a tr ) < 1. Figure 5 shows this expression as a function of atr/oWn- Taking 
a m in = O.OlxlO -6 m and a tr = 4.5xl0~ 14 v 2 m with v in km/s, we get 

. 4.0xl0 7 9.7xl0 10 

<W - "^73 — ( 38 ) 

Using the values from above, namely, A = 8.9xl0~ 14 v 2 m from equation (17), n = 4.0xl0 5 m~ 3 from equa- 
tion (8), C = 2.0xl0" 15 m n / 3 , and R = 1738xl0 3 m, but now taking a min = O.OlxlO^m for the smallest dust 
grain radius considered here, we get 

= 6.4xl0- 15 v 4 ■ C{v), v > 474 km/s (39) 

the difference between this and equation (30) being the presence of the factor C(v). Thus, we have for the UV 
model: 

rft 

v < 474 km/s = v cr i t 
», v > 474 km/s = v crlt (40) 

Figure 6 shows the ratio of neutral solar wind flux to solar wind flux from equation (40) on the left y-axis 
and solid line assuming the minimum dust grain radius is 0.01/zm. The right y-axis and dotted line shows the 
corresponding count rate assuming the same parameters as in Figure 4. At the highest solar wind speeds, a 
LENA-like instrument would observe a count rate of over 100 s _1 . For all except the very slowest solar wind 
speeds, such an instrument would register a count rate > 1 s _1 . 

Note that, although the solar wind is typically about 96% protons by number, smaller neutral atom fluxes 
of minor species, notably helium and oxygen, would also be created by the solar wind interaction with lunar 
dust. 

The validity of the approximation z = x 2 /2R used in equation (20) is dependent on the maximum dust 
scale height in comparison to the lunar radius. This scale height, C/c?J^ nl is only 9.3 km for the "visible" light- 
scattering dust model discussed in the previous section corresponding to a min = 10~ 7 m and so the approximation 
is valid. However, for the "UV" scattering model, this scale height is 4300 km (10 8 / 3 = 464 times larger) which is 
about 2.5 times larger than the lunar radius of 1738 km. However, the overall accuracy of the approximation will 
be bounded by how well the spatial integral using z — x 2 /2R compares to the integral, evaluated numerically, 
without this approximation. The spatial integral evaluated numerically as compared with the approximation in 
equation (25) given by 7.4xl0~ 5 a~ 4 / 3 as a function of dust radius shows that the maximum discrepancy occurs 
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at the smallest dust radius, as expected, and is only about 20%. Worth noting is that the approximation actually 
underestimates the spatial integral such that the actual flux of lunar dust-generated neutral solar wind will be 
larger than predicted by the analytic solution. 

7. Discussion 

It has been the purpose of this paper to examine neutral solar wind created by lunar exospheric dust. 
However, as mentioned in section 1, neutral solar wind also results from solar wind charge exchange (SWCX) 
with exospheres [e.g. Collier ct al., 2005] and this process also occurs at the Moon. If we assume an argon 
atmosphere with a surface density of uat = 10 5 cm -3 , a scale height of Hat = 40 km and a canonical collision 
cross section gat = 10~ 15 cm 2 [Stern, 1999], one can estimate a ratio of neutral solar wind to solar wind flux 
viewing normal to the surface ^nsw/^sw ~ nArOArHAr = 4xl0~ 4 . Assuming limb viewing, ^nsw/'I'sw will be 
~ 10~ 3 . Referring to Figure 6, at higher solar wind speeds the lunar dust-generated neutral solar wind flux will 
likely exceed that due to charge exchange while at slower solar wind speeds the opposite will be the case. 

However, even at slower solar wind speeds there are a number of factors which would allow a separation 
between charge-exchange and lunar dust-generated neutral solar wind: (i) Dust-generated energetic neutral atoms 
(ENAs) would not be as beam-like as the charge- exchange generated ENAs because penetrating dust grains will 
scatter the resulting neutrals in angle; (ii) Lunar dust-generated neutral solar wind would contain species other 
than just hydrogen because penetrating dust can neutralize high charge state ions; (iii) Lunar dust-generated 
neutral solar wind would have a significantly wider energy distribution than charge- exchange generated neutral 
solar wind because the process of penetrating dust grains involves energy loss and straggling; (iv) Finally, the 
relative flux of dust-generated neutral solar wind increases dramatically with solar wind speed compared with 
charge exchange with neutral exospheric atoms due to the solar wind speed-dependent cross-section [e.g. Rucinski 
et al., 1996]. 

Along similar lines, the ratio of the neutral solar wind to ionized solar wind due to charge transfer in the 
interplanetary medium (i.e. without an exosphere) is expected to be 10 -5 — 10~ 4 [e.g. Holzer, 1977; Gruntman, 
1994], so that the lunar dust-generated neutral solar wind signal is comparable to this population at slow solar 
wind speeds but dominates it at higher solar wind speeds. Furthermore, because the lunar dust-generated neutral 
solar wind is anticipated to be a very local phenomenon, changes in its flux will be more highly correlated than 
the "background" neutral solar wind with ionized solar wind observations from an upstream monitor within a 
few tens of Earth radii from the location of the Moon [e.g. Collier et al., 1998; Richardson and Paularena, 2001]. 

There were a number of simplifying assumptions that went into the analytic results in the previous sections. 
Among these is the use of data on the penetration of protons through silicon which were readily available due to 
the commercial importance of implanted silicon. However, a more accurate approach might be to use data for 
the penetration of protons through silica [Eder et al., 1997]. Although the data are not as complete as those for 
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protons through silicon, the data suggest the range is higher for protons through silica. We also assumed that 
lunar dust is spherical when in reality its shape is more complex. These assumptions may lead to a neutral solar 
wind flux underestimated by a factor of four or more. 

Some of the other assumptions are that protons execute straight-line trajectories through a dust grain. 
Although this is a reasonable first approximation, some protons will back-scatter from the dust grain while those 
that penetrate will both scatter in angle and lose energy. Furthermore, the penetrating protons will sputter 
both ions and neutral atoms from the dust grains. None of these effects are well-understood. For example, the 
sputtering yield formulae may not be accurate for small spherical dust grains because they are obtained from 
experiments on flat solids. Furthermore, transmission sputtering (sputtering on exit) may be important, but 
has not been investigated in the regime of interest here [e.g. Dwek and Arendt, 1992]. However, the sputtered 
products will have relatively low energies (~10 eV) in comparison to the solar wind [Johnson and Sittler, 1990; 
Elphic ct al., 1991]. Consequently, they should be easily separated out from the more energetic neutral solar 
wind component. 

The energy distribution will also be affected by the solar wind-dust interaction. Ions traversing through 
matter lose energy by a sequence of stochastic processes which result in an energy loss and dispersion. Typically 
at incident energies around 4 keV, proton beams traversing solid foils between 0.01 and 0.03 ^m thickness will 
emerge with ~3 keV with energy dispersions of the order of a few hundred eV [Figueroa et al., 2007]. So, the 
observed energy distribution of the lunar dust-generated neutral solar wind will have lower energies than that of 
the solar wind and will be substantially hotter. Also note that energy arguments suggest that a 1 keV proton 
depositing all its energy into a 0.01 /xm dust grain will not destroy the dust grain or significantly affect the overall 
integrity of the grain. 

8. Conclusions 

Dust is a recognized planetary hazard which will have an adverse impact on robotic and human exploration 
of the Moon and Mars. In this paper, we showed that based on our current understanding of lunar exosphcric 
dust characteristics at the terminator, a significant and easily observable flux of energetic neutral solar wind 
hydrogen, at times higher than 10~ 3 the solar wind flux even when conservatively estimated, will be present 
in the terminator region. Thus, neutral atom observations would provide information on the exospheric dust 
distribution at the Moon. For example, the distribution of small grains, which are easier for the solar wind 
to penetrate, could be determined if one assumes the altitude profile follows the Murphy and Vondrak [1993] 
distribution used in the paper. Furthermore, other model distributions could be tested [e.g. Stubbs et al., 2007d]. 
Consequently, neutral atom observations would be highly complementary to optical techniques for characterizing 
lunar dust and should be considered an important component of any lunar robotic mission designed to assess 
the characteristics of this hazard. 
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It has been speculated that activities associated with in-situ resource utilization (ISRU) could increase the 
amount of exospheric lunar dust by several orders of magnitude [Vondrak, 2007]. Given that these activities 
may very well occur at the polar regions where the terminator is omnipresent and that the ratio of neutral solar 
wind flux to solar wind flux increases linearly with dust concentration (i.e. no), such activities may very well 
neutralize a large fraction of the solar wind. 

Retarded space weathering due to magnetic anomalies shielding the lunar surface from the solar wind has 
been proposed as an explanation for high albedo markings [e.g. Hood and Williams, 1989]. Because neutral 
solar wind will not be deflected by magnetic anomalies, a substantial flux of neutral solar wind could alter space 
weathering processes. 
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Figure Captions 

Figure 1. The dust concentration in m~ 3 as a function of altitude above the lunar surface at the terminator 
for dust grain radii from 0.1-6 /im (solid line) from Stubbs et al. [2007d] (based on a numerical integration of 
the Murphy and Vondrak [1993] model) and the approximation (dashed line) given by equation (16). 

Figure 2. The range of hydrogen through silicon as a function of hydrogen speed over the typical range of 
solar wind speeds. Note that even at the highest solar wind speeds, ^1000 km/s, solar wind protons will not all 
penetrate spherical dust grains with radius 0.1 /im. However, at solar wind speeds slightly greater than average 
(^440 km/s), solar wind protons will completely penetrate 0.01 /zm grains. 

Figure 3. Geometry of neutral solar wind viewing at the lunar terminator. The Sun is to the left and the 
solar wind flows antisunward with speed v. The lunar radius, 1738 km, is designated R. A point x along the 
line-of-sight corresponds to an altitude z above the lunar surface. 

Figure 4. Prediction for the ratio of neutral solar wind flux to ionized solar wind flux considering only the lunar 
dust population at the terminator with grain radii in excess of 0.1 /xm (left axis, solid line). Also, the neutral 
solar wind count rate that would be observed by an instrument like IMAGE/LENA assuming a solar wind flux 
of 3xl0 8 /cm 2 /s and an efficiency of 2xl0~ 4 as a function of solar wind speed (right axis, dashed line). 

Figure 5. The correction factor C(at r ) as a function of a tr /a m i n . C(at r ) modifies the neutral solar wind flux 
calculation to take into account the effect of the solar wind fully penetrating the dust grain. Recall that a tr is a 
function of solar wind speed. 

Figure 6. Prediction for the ratio of neutral solar wind flux to ionized solar wind flux considering the lunar 
dust population at the terminator with grain radii in excess of 0.01 /zm (left axis, solid line). Also, the neutral 
solar wind count rate that would be observed by an instrument like IMAGE /LENA instrument assuming a solar 
wind flux of 3xl0 8 /cm 2 /s and an efficiency of 2xl0~ 4 as a function of solar wind speed (right axis, dashed line). 
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